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ABSTRACT 


The  parametric  adaptive  matched  filter  (PAMF)  for  space-time  adaptive  processing 
(STAP)  has  been  an  interesting  research  area  in  the  airborne  radar  data  processing 
community  over  the  last  decade.  Starting  from  providing  a  complete  formulation  of 
PAM  F,  this  report  assesses  the  performance  of  PAM  F  using  simulated  data  as  well  as 
real  airborne  data  collected  by  the  Multi-channel  Airborne  Radar  Measurements 
(MCARM)  system  with  use  of  results  from  the  conventional  STAP  as  a  benchmark.  A 
modified  PAMF  approach  using  both  forward  and  backward  predictions  is  proposed 
to  eliminate  the  dimensionality  loss  of  the  method.  Finally  the  operational  counts  and 
computati  onal  savi  ngs  of  PA  M  F  are  al so  esti  mated . 
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Parametric  Adaptive  Matched  Filter  and  itsModified 

Version 


EXECUTIVE  SUMMARY 

How  to  quickly  process  phased  array  radar  data  for  small  target  detection  is  a  key 
issue  of  air  surveillance  radar  systems.  Much  research  effort,  in  particular  that 
sponsored  by  the  Air  Force  Research  Laboratory  (AFRL),  has  been  put  into  developing 
the  parametric  adaptive  matched  filter  (PAMF)  for  space-time  adaptive  processing 
(STAP)  over  the  last  decade.  This  report  traces,  examines  and  evaluates  this  technique 
in  great  detail.  The  contribution  of  this  report  is  manyfold.  First  a  full  mathematical 
formulation  of  PAMF  is  provided.  Previously  published  papers  on  this  topic  often  omit 
details  of  the  mathematics,  possibly  due  to  space  limitation,  making  it  difficult  for 
readers  to  follow  through.  This  report  provides  an  independent  derivation  for  the 
PAMF  algorithm.  With  the  details  of  the  mathematics  provided,  the  concept, 
advantages  and  shortcomings  of  the  filter  should  beunderstood  more  thoroughly. 

Secondly  the  performance  of  PA  M  F  i  s  ful  ly  assessed  by  applyi  ng  the  al  gorithm  to  three 
different  types  of  radar  data  including  the  data  generated  using  a  generic  radar  model, 
data  generated  using  high  fidelity  airborne  radar  simulation  software,  RLSTAP  and 
real  airborne  radar  data  collected  by  the  M  ulti-Channel  Airborne  Radar  Measurements 
(M  CARM )  system.  The  results  are  compared  to  benchmarks,  the  results  of  space-time 
adaptive  processing  (STAP). 

Other  authors,  eg,  Roman  et  al  (2000)  have  observed  that  PAMF  is  robust  in  the 
situation  of  having  reduced  sample  data  support,  but  the  reason  remained 
unexplained.  This  paper  mathematically  explains  the  reason.  There  is  an  equivalent 
covariance  matrix  for  PAMF,  whose  dimension  is  much  smaller  than  the  covariance 
matrix  for  the  conventional  STAP.  The  number  of  the  dominant  eigenvalues  of  the 
equivalent  covariance  matrix  for  PAMF  is  also  much  smaller  than  that  for  the 
conventional  STAP,  so  that  less  sample  data  are  required  for  estimating  the  equivalent 
covariance  matrix.  From  the  processing  point  of  view  PAMF  uses  both  fast-time  and 
slow -time  averaging  processing  in  forming  the  equivalent  covariance  matrix,  whereas 
the  conventional  STAP  only  utilises  fast-time  averaging. 
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Another  contribution  is  the  estimation  of  operational  counts  (ops)  for  the  PAMF 
algorithm.  Computational  savings  of  PAMF  in  comparison  to  the  conventional  STAP  is 
also  provided. 

Finally  a  modified  PAMF  approach  is  proposed.  The  existing  PAMF  approach  only 
uses  forward  prediction  parameters,  leading  to  a  dimensionality  loss,  which  in  turn 
reduces  the  Doppler  resolution  and  coherent  processing  gain.  In  the  case  where  the 
order  of  the  filter  is  high  and  the  number  of  pulses  in  a  coherent  process  interval  (CPI) 
is  small,  the  loss  in  coherent  processing  gain  and  the  loss  in  Doppler  resolution  could 
be  significant.  The  proposed  PAMF  approach  uses  both  forward  and  backward 
prediction  parameters  to  eliminate  the  dimensionality  loss. 

Based  on  this  study,  the  PAMF  algorithm  has  two  advantages  over  the  conventional 
STAP  algorithm: 

•  The  requirement  of  operational  counts  (ops)  is  reduced.  Computational 
savings,  compared  to  the  STAP  algorithm,  vary  from  less  than  20% for  a  system 
with  a  small  number  of  pulses,  a  small  number  of  array  elements  and  a  large 
number  of  range  bins,  to  more  than  90%  for  a  system  with  a  large  number  of 
pulses,  a  large  number  of  array  elements  and  asmall  number  of  range  bins. 

•  The  PAMF  algorithm  is  much  more  robust  than  the  STAP  algorithm  in  the  case 
where  the  set  of  sample  data  used  for  generating  the  covariance  matrix  is 
relatively  small. 

Problems  of  the  PAMF  process  include 

•  Because  of  the  dimensionality  reduction,  the  equivalent  temporal  independent 
measurements  in  a  CPI  are  reduced,  which  in  turn  reduces  the  Doppler 
resolution  and  coherent  processing  gain.  The  proposed  PAMF,  which  uses  both 
forward  and  backward  prediction  parameters,  solves  the  dimensionality 
reduction  problem  and  maintainstheDoppler  resolution,  and  processing  gain. 

•  The  cl  utter  profile  against  range  bin  using  the  PAMF  processing  does  not  seem 
to  follow  the  pattern  of  the  original  clutter  profile.  Whether  this  imposes  a 
problem  for  target  detection  warrants  further  investigation. 

All  simulations  have  suggested  that  using  a  sliding  window  to  sample  the  covariance 
matrix  to  process  each  range  bin  does  not  statistically  improve  results.  Therefore  a 
preferable  option  isto  useafixed  window  to  samplethe  covariance  matrix  and  process 
a  block  of  range  bins,  which  significantly  saves  computational  cost. 

Overall  the  amount  of  ops  required  for  the  PAMF  algorithm  is  still  huge  considering 
the  time  frame  available  for  airborne  radars.  Flow  to  further  reduce  ops  demand  in 
processing  phased  array  data  warrants  future  research. 
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1.  Introduction 

In  a  scenario  using  phased  array  radar  for  detecting  moving  targets  whose  signals  are 
embedded  in  undesired  signals  including  clutter  (echoes  from  the  Earth's  surface), 
broadband  jamming  and  thermal  noise,  space-time  adaptive  processing  (STAP)  yields 
an  optimal  signal-to-interference-and-noise  ratio  (SINR).  However,  there  are  two 
critical  issues  in  real-time  implementation  of  STAP  in  airborne  radar  systems.  First  for 
a  range  cell  under  test  (CUT),  data  from  a  large  number  of  adjacent  range  cells  are 
required  in  order  to  adaptively  compute  a  reliable  covariance  matrix  of  undesired 
signals.  The  available  sample  cells  are  sometimes  fewer  than  the  required  for  accurate 
estimation  of  the  covariance  matrix.  This,  however,  may  be  solved  by  assuming  the 
structure  of  the  covariance  matrix  to  be  known  (Steiner  and  Gerlach,  1998,  Steiner  and 
Gerlach,  2000,  Gerlach  and  Picciolo,  2003,  Bresler,  1988).  The  second  and  more  critical 
issue  is  the  computational  requirement,  which  in  fact  becomes  a  processing  bottleneck 
at  current  computer  speeds.  Numerous  algorithms,  aimed  at  reducing  the 
dimensionality  of  the  covariance  matrix  in  both  the  spatial  and  temporal  domains, 
have  been  proposed  to  minimise  the  computational  demand  to  satisfy  the  real-time 
requi  rement  while  mai  ntai  ni  ng  the  coherent  processi  ng  gai  n  at  or  cl ose  to  the  opti  mum 
level  (reduced  dimensionality  also  requires  fewer  sample  cells)  (Wang  and  Cai,  1994, 
Ward  and  Kogon,  2004,  Wang  et  al,  2003). 

M  uch  research  effort,  notably  that  carried  out  and  sponsored  by  the  Air  Force  Research 
Laboratory  (AFRL),  has  been  done  towards  developing  and  evaluating  the 
performance  of  the  parametric  adaptive  matched  filtering  (PAMF)  approach  in  recent 
years.  Many  papers  have  been  published  on  this  topic  in  the  last  decade  (Roman  et  al, 
2000,  Michels  et  al,  2003,  Rangaswamy  and  Michels,  1997,  Rangaswamy  et  al,  1995, 
Robey  et  al,  1992,  Roman  et  al,  1997).  The  advantages  of  PAMF  include  a  reduced 
sample  requirement  in  estimating  parameters  and  a  reduced  computational 
requirement  (Roman  etal,  2000). 

Building  from  a  mathematical  derivation  of  the  PAM  F  (Sections  2  to  4)  that  is  missing 
from  the  I  iterature  (space  I i  mitati on)  the  contri  buti ons  of  this  report  are: 

a)  An  assessment  of  the  performance  of  the  PAMF  compared  to  the  traditional 
matched  filter  (eg,  STAP)  in  Section  5; 

b)  An  explanation  of  the  observed  robustness  of  the  PAMF  when  applied  to 
reduced  data  samples; 

c)  Estimations  of  the  computational  load  of  the  PAMF  approach  in  Section  6; 

d)  In  Section  7  a  modified  PAMF  approach  that  eliminates  the  dimensionality  loss 
of  the  PAM  F  to  recover  the  processing  gain  and  Doppler  resolution  lost  by  the 
PAMF. 
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2.  Matched  Filtering  Processing 

The  mathematical  notation  used  in  this  report  follows  the  convention:  boldface 
uppercase  and  lowercase  letters  represent  matrices  and  vectors,  respectively, 
superscripts  T  and  H  denote  transpose  and  Hermiti  an  transpose,  respectively,  and  C 
and  (R  stand  for  the  complex  and  real  number  fields,  respectively.  For  instance, 
xe€“xi  denotes  x  to  be  an  MN -element  complex-valued  column  vector,  and 
AeCWxW  denotes  A  to  be  an  MN x  MN  -el ement  complex-valued  matrix.  The 
symbol  0  denotes  the  Kronecker  matrix  product,  and  E{-}  expectation. 

Let  a  snapshot  {x(n)|«  =  0,l,---,M-l}  with  x(n)e€m  be  the  complex  valued 

sequence  at  the  output  of  a  linear  equispaced  N  -element  array  after  demodulation  and 
pulse  compression  sampling,  corresponding  to  the  return  from  a  single  range  cel  I  over 
an  M  -pulseCPI.  If  undesired  signals  xu  are  Gaussian  distributed,  itiswell  known  that 
the  optimal  MF  weighting  vector  is  (Compton,  Jr,  1988  and  Ward,  1994), 

WMF=rK1e  (i) 

where  y  is  an  arbitrary  scalar,  which  without  loss  of  generality,  we  will  let  equal  1. 
R1(  =  E{xux" }  e  cMNxMN  is  the  covariance  matrix  of  undesired  signals.  The  undesired 
signals  xu  may  consist  of  clutter,  thermal  noise  and  possibly  noise  jamming.  The 

vector  e  e  CM7Vxl  is  the  desired  signal  spatial -temporal  steering  vector,  whose  explicit 
expression  may  be  written  as, 

e  =  e(/,)®e(/J  (2) 

where  ftd  and  fts  are  the  normalised  target  Doppler  and  spatial  (azimuth) 
frequencies,  respectively.  Vectors  e(ftd )  and  e(/te)  are  the  desired  target  temporal 
and  spatial  steering  vectors,  respectively,  defined  by, 

e(A/)  =  ^=[l  exp(y2^)  •••  exp(j2x(M  -  \)ftd)]r  (3) 

and 

e(/r,)  =  ^=t1  expO‘2^,)  •••  exp(j2^-(A-l)/te)f  (4) 

Using  block  component  notation,  (2)  may  be  written  as, 
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n  =  -1  (5) 

The  output  of  the  MF  processor  is  the  product  of  the  Fiermitian  transpose  of  the 
opti  mal  wei  ghti  ng  vector  ti  mes  the  data  snapshot, 

y  =  wmfx  (6) 

The  coherent  processing  gain  of  the  MF  processor  approaches  the  upper  limit, 
101og10(MAT)dB,  for  deterministic  signals  whose  bearings  differ  from  the  jamming 
bearings,  and  whose  Doppler  frequencies  differ  from  the  clutter  Doppler  frequency  in 
the  radar  look  direction. 


Given  a  target  signal  ere ,  the  SI  NR  is  defined  as  the  ratio  of  the  output  target  power  to 
the  output  i  nterference  and  noise  power, 


SINRmf 
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=  £e"R,> 


(7) 


where  £  =E{\  a  |2} . 

The  response  of  the  weighting  vector  to  the  azimuth  and  Doppler  is  referred  to  as  the 
adapted  pattern  and  defined  as  (Ward,  1994) 

PMF=\<AftsJld)\2  (8) 


where  e(fts,fld)  is  a  general  steering  vector  with  Doppler  frequency  ftd  varying  from 
-0.5  to  0.5,  and  spati  al  frequency  fts  (cosi  ne  of  the  azi  muth)  varyi  ng  from  -1  to  1. 


For  an  unknown  signal  amplitude,  a  detection  test  statistic  (DTS)  for  the  MF  process 
has  been  proposed  as  the  ratio  of  the  output  power  to  the  power  of  interference  and 
noise  (Robey  et  al,  1992), 


MF 


I  w  H  X  I 
I  ™  mfa  I 


w  R  w 

MF  it  ™  MF 


e^R^x  | 
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(9) 


Comparing  (7)  and  (9)  we  can  seen  that  the  MF  DTS  is  actually  the  SI  NR  of  the  signal 
x .  In  operations,  when  XMF  exceeds  a  threshold,  target  presence  is  declared. 
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The  DTS  AMF  may  be  expressed  in  an  alternative  way  by  utilising  the  LDU  block 
decomposition  of  R;( .  Because  R1(  is  a  positive  definite  and  Hermitian  matrix,  it  can 
be  decomposed  in  the  form  of  RJ(  =  ADA22 .  In  this  form,  \e£MNxMN  \s  a  lower 
block-triangular  matrix  with  M  identity  matrices  lN  along  its  main  block  diagonal, 
and  D e £MNxMN  \s  a  block-diagonal  matrix  with  matrices  D0i)eC", 
n  =  along  its  main  block  diagonal.  Substituting  R~‘  =  A_//D“1/2D“1/2A_1 

in  (9)  results  in, 


|  (D~1/2A~1e)//(D~1/2A~1x)  |2  _  |  (D~1/2u)//(D~1/2£)  |2 
(D-1/2A-1e)//(D^1/2A“1e)  ~  (D_1/2u)^(D“I/2u) 
sgv|2  |  |2 

sHs  Zf=o1s//(n)s(n) 


(10) 


where  u  =  A  e,  s  =  D  1  2u  ,  e  =  A  lx  and  v  =  D  1/2e  .  Let  A  1  bedenoted  as, 


Ijv  0  •••  0 

Af(l)  lN  -  0 

Am_i(M—  1)  Am_,(M  — 2)  •••  Ijy 


(ID 


The  block  elements  s(n)  and  \(n)  can  be  written  as, 


e(n)  =  x(n)+Z^(k)x(n-k)  n  =  -1  (12) 

k= i 

v(«)  =  D^1/2  («)£(«)  n  =  -1  (13) 


It  can  be  proven  that  with  or  without  target  signal  presence,  v  has  identity  covariance 

matrix.  Therefore,  A-1  implements  block  whitening  and  D~1/2  implements  block 
unity-variance  normalisation.  After  the  whitening,  vectors  in  the  sequence 
{je(w)| «  =  0,l,--,M-l}  are  uncorrelated  in  pairs  and  the  covariance  of  s(n)  is  D (n) . 
After  further  block  unity-variance  normalisation,  vectors  in  the  sequence 
\y(n)\n  =  -l}  are  uncorrelated  in  pairs  and  each  has  identity  covariance.  In 

another  words,  the  undesired  signal  xu  that  originally  satisfies  the  complex  normal 
distribution,  with  zero  mean  and  covariance  R„ ,  i.e.,  xu~CN( 0,  R„),  is  now 
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transformed  into  v  that  satisfies  the  complex  normal  distribution,  with  zero  mean  and 
identity  covariance,  i.e.,  v  ~  CN( 0, 1) . 

Equation  (12)  may  be  viewed  as  a  result  of  a  moving-average  (M  A)  filter  in  block  mode 
(Roman  et  al,  2000).  The  signal  x(n) ,  n  =  0,1,---,M-1,  is  estimated  by  the  forward 

linear  prediction  using  x(n-k),  k  =  1, •••,«,  with  the  parameters  A f(£),  k=l,--,n. 
The  prediction  error  is  z(n ) . 


3.  Parametric  Matched  Filter 

This  Section  formulates  the  parametric  matched  filter  (PMF).  Roman  et  al  (2000) 
mention  that  some  simplifications  in  theMF  DTS,  AMF ,  lead  to  the  PMF  DTS,  APMF ,  to 
be  given  below.  In  fact  there  is  no  direct  mathematical  connections  between  AFM  and 
Apmf  ,  that  is,  we  cannot  derive  APMF  from  AMF ,  although  we  may  express  both  AMF 
and  Apmf  in  similar  forms. 

Assume  the  data  sequence  {x(n)\n  =  p,p  +  -l}  is  estimated  using  the  pth  order 

forward  autoregressive  (AR)  parameters  A Hf(k),  k  =  l,--,p.  In  this  process  the 
forward  linear  prediction  x(n)  for  sample  x(n)  is  a  linear  combination  of  x(n-k), 
&  =  1,- as  (Marple,  Jr,  1987), 


x(ri)  =  -YJAHf(k)x(n-k)  n  =  p,p +  (14) 

k= 1 

Theforward  linear  prediction  residual  is, 

z(n)  =  x(n)-x(n)  =  X  AHf(k)x(n-k)  n  =  p,p  +  \,---,M -1  (15) 

k= 0 

where  A^(0)  =  Iw.  The  count  n  of  e(n)  in  (15)  starts  from  p ,  which  is  not 
convenient.  By  shifting  the  count  n  of  e(n)  to  start  from  zero,  (15)  can  be  rewritten  as, 

£(«)  =  Z  AHf  (k)x{n  -  k  +  p)  n  =  0,l,---,M  -  p-l  (16) 

k= 0 

Although  (16)  and  (12)  use  similar  symbols,  this  is  adopted  to  avoid  using  too  many 
symbols.  Note  that  the  subscript  /  of  A Hf(k)  in  (16)  simply  means  the  forward 
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prediction.  Unlike  in  (12),  A^(A:)is  now  independent  of  n .  Also  A ^(k)  in  (16)  in 
general  has  nothing  mathematically  to  do  with  Af  (k)  in  (12).  This  will  become  clear 
after  the  derivation  of  A ^(k)  is  given  in  thefol lowing  subsection. 

Similarly  we  may  define  the  steering  residual  ueC(M_p)7Vxl  as, 

u(n)=ZAHf(k)e(n-k  +  p)  n  (17) 

k= 0 

It  is  worth  noting  that  the  dimensionality  of  the  pth  order  £  and  u  is  reduced  to 

^(M-p)Nxl 

Equations  (16)  and  (17)  can  be  written  in  more  compact  forms  as, 

£(0) 

£(1) 

£  (M  -  p  - 1) 

ii(0) 

u(l) 

u(M  -  p  -1) 

where  the  block  matrix  B  e  <^m~p)n*mn  has  the  form  of 

A  Hf(p)  A*(p-l)  -  lN  0  -  0 

0  A  Hf{p)  A  ^  (p  - 1)  -  lN  0  -  0 

0  -  0  A  Hf{p)  Aj  (p  - 1)  -  lN 

Let  R e(n)  =  E{£(n)£H(n)}  denote  the  temporal  residual  covariance  matrix  which  is 
Hermitian  and  positive  definite.  TheLDU  decomposition  of  R  E(n)  is, 

R»  =  L»D»L"(/i)  (21) 
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where  L e(n)  e  CNxN  is  a  lower  triangular  matrix  and  its  main  diagonal  elements  are 

unity  and  DE(n)  <eUNxN  \s  a  diagonal  matrix  with  its  all  diagonal  elements  greater 
than  zero.  The  transfer  matrix  T (n)  functioning  as  both  the  whitening  and  unity- 
variance  normalisation  for  s(n)  isthen, 


T(n)  =  D;1/2(n)L->) 


Finally  the  PM  F  DTScan  be  defined  as, 


A 


PMF 


|s*v|: 

“77 


Zf=/-,s//(n)s(n) 


(22) 


(23) 


where  s(«)  =  T(«)u(«)  and  v(n)  =  T(n)e(n) .  Equation  (23)  can  be  expressed  in  a 
compactform  with  theuseof  input  signal  x  and  steering  vector  e  as, 


A 


PMF 


|  s^v  |2 

sHs 


|  e^CpX  |2 

eHC  e 


where 


cP  =b"r;1b 


(24) 


(25) 


R 


£ 


R,(0) 

R,(i) 


KE{M-p-\) 


(26) 


Therefore,  analogous  to  STAP,  we  may  define  the  PMF  weighting  vector 


spMNxl 

w  pmf  e  as, 


w 


PAMF 


=  rCPe 


(27) 


The  corresponding  SIN  R  and  adaptive  pattern,  respectively,  are, 


S/NR 

and 


4t  I  W  PMFe  I 


. H , 


PMF  \yH  R  w  e"C  R  C  e 

w  PMF^u  ”  PMF  e 


=  &HCpt 


(28) 
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PpMF  =\<MF*{ftsJtd)\2  (29) 

Comparing  (24)  and  (9)  or  (28)  and  (7),  we  see  that  the  matrix  CP  in  (24)  and  (28) 
functions  as  the  matrix  in  (9)  and  (7). 


3.1  Determination  of  Parameters  in  Parametric  M  atched  Filtering 

In  the  previous  subsection  we  have  derived  for  the  PMF  algorithm  all  formulae  except 
the  filter  parameters  AHf(k),  k  =  l,---,p.  Different  methodologies  lead  to  different 

algorithms  for  finding  the  parameters  (Marple,  Jr,  1987).  Here  we  follow  the  least 
squares  method,  i.e.,  the  covariance  method  (Marple,  Jr,  1987).  Roman  et  al  (2000)  has 
reported  that  in  the  context  of  radar  signal  processing,  the  least  squares  method 
usually  outperforms  other  methods,  such  astheStrand-Nuttall  method.  The  derivation 
bel  ow  i  s  i  ndependent  work  of  the  author. 

The  goal  of  the  least  squares  method  is  to  find  parameters  A  f(k),  k  =  \,---,p,  which 

produce  a  minimum  squared  residual,  ehe.  To  find  the  minimum  squared  residual, 

8 

we  need  to  solve  the  simultaneous  equations  — - — ehe  =  0,  k  =  \,--,p  (note  that 

8AHf  (k) 

A^(0)  =  IN  is  constant).  The  derivative  of  a  scalar  function  / (Q)  e  <Clxl  with  respect 
to  Q  e  <C"!X'!  is  an  m  x  n  matrix  (Van  Trees,  2002,  p.  1399), 


a 

8Q 


/( Q)  = 


df  /8qn 
df/dqml 


df  <dqln 

df!8qmn 


(30) 


where  qik,  i  =  !,■■■, m,  k  =  \,--,n,  is  the  element  of  Q.  After  tedious  but 

d 

straightforward  algebra,  the  equation  satisfying  — — — ehe  =  0,  k  =  1,  •••,/>,  is, 

3 A.  j-  (k') 


K  *  R/i. 


aap ) 


ATI) 


LN 


=  0eC 


pNxN 


Therefore  we  find  the  parameters  as, 


(31) 
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A  f(p) 

-R^Ryj 

Ay-  =  Ay(l)  =  (32) 

_  Ijv 

_Ay(0)_ 

where  Ry  eC^  and  Rfl  eCpNxN  are  sub-block  matrices  of  Ry  e  <c<^+i)zvx(P+i)zv 
given  by, 

R^r  R/i 

Rf=  — - -  =E{XfX Hf)  (33) 

Rio  Rn_ 

x(0)  x(l)  •••  x(M-p- 1) 

X/=  X<1)  X(2)  X(M-p)  (34) 

_x(p)  x(/?  +  l)  •••  x(M-l) 

If  the  covariance  matrix  R„  is  known,  R;  can  be  formed  by  summing  up  M-p 
block  matrices  of  RJ(  usingamoving  (p  +  \)N  x  (p  + 1)  N  window  diagonally  from  top 
to  bottom,  which  may  be  expressed  as, 

M-p- 1 

R/-  =  Z  R u(k:k  +  p)  (35) 

k= 0 

where  R  u(k:k  + p)eC{p+l)Nx(p+l)N  is  a  block  sub-matrix  of  R;( ,  with 
Rt,  ( k,k )  e  CiVxA  as  its  top  left  block  element  and  R£,  (k  +  p,k  +  p)e  CNxN  as  its  bottom 
right  block  element.  Figure  1  depicts  a  situation  where  Ry  is  formed  by  adding  four 
sub-block  matrices  of  R„ ,  corresponding  to  the  case  of  p  =  2  and  M  =6.  It  can  be  seen 
that  the  PM  F  algorithm  does  not  use  all  block  elements  in  R„  to  form  R  f .  This  may 
be  considered  as  an  information  loss  since  the  block  element  R„(z',&)  stands  for  the 

correlation  between  data  collected  by  the  z'th  and  kt h  pulses  and  its  neglect  means  that 
the  correl  ati  on  i  s  treated  as  zero. 
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the  case  of  p  =  2  and  M  =  6.  Some  of  the  block  elements  of  Ru  are  not  used  In  forming  R;. 

4.  Parametric  Adaptive  M  atched  Filter  (PAM  F) 

In  general  the  covariance  matrix  of  undesired  signals  Rt,  is  unknown,  so  it  has  to  be 
estimated  adaptively  using  range  samples.  The  often  used  method  is  the  diagonally- 
loaded  sample  matrix  (DL-SM)  method  (Carlson,  1988,  Ward  and  Kogon,  2004),  as, 

R»=^SV?+<5I  (36) 

A  k= 1 

where  the  subscri  pt  k  denotes  range  bins.  8  is  a  small  value  of  the  order  of  the  system 
thermal  noise  level. 

The  transfer  matrix  T(«)  n  =  0,1, p- 1  given  by  (22)  is  time-dependent.  Roman 

et  al  (2000)  suggest  that  a  single  time-independent  transfer  matrix  T  may  replace  the 
time-dependent  transfer  matrix  T(n)  in  the  adaptive  process.  The  time-independent 

transfer  matrix  T  is  obtained  through  a  ti me-averagi ng  process  as, 

—  1  M~p~l  1  „ 

R /,=— -  I  R»  =  — - XHfRfXf  (37) 

M  -  p  n= o  M  -  p 

(38) 
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t  =  d;1/2l-; 

ThePAMF  DTS  is  defined  as, 


A 


_  |  s^v  | 

PAMF  ~  h 

S  S 


lgj>)vW  I 

I%rx*H(n)  s(») 


(39) 


(40) 


Expressions  (40)  and  (23)  are  the  same,  but  the  difference  between  APMF  and  Apamf  is 
that  Ru  and  T (n)  n  =  - p-\ ,  are  used  for  PMF  whereas  Rt,  and  T  are 

used  for  PAMF. 

Equation  (40)  can  be  written  in  a  more  compact  matrix  notation  with  input  signal  x 
and  steeri  ng  vector  e  (Michelsetal,  2003).  In  order  to  do  so,  we  write, 


v(0) 

T 

£(0) 

v(l) 

= 

T 

£(1) 

v(M  -  p-l) 

T 

e(M  - p- 1) 

This  equation  can  be  written  as, 

v  =  {lM_p  0  T)e  =  {lM_p  0  T)Bx  (42) 

Similarly, 

s  =  (IM_p0T)u  =  (I M_p  0  T)Be  (43) 

Therefore, 

swv=eHC?j4x  (44) 

sffs  =  eHCP/  (45) 

where 

CPA  =Vh{Im_p  ®TH)(IM_p  0 T)B 

=  B//(Im_p0T//T)B  (46) 
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I  nserti  ng  (44)  and  (45)  i  n  (40),  we  can  express  the  PA  M  F  DTS  as, 


A 


PAMF  ~ 


PHC  Yl2 

eHC  e 


(47) 


Similarly  the  matrix  CPA  in  (47)  functions  as  R,,1  in  (9).  The  PAMF  weighting  vector 
w Pamf  e  CA/X  defined  implicitly  in  (47)  is, 


w pamf  ~  yCpAe 

The  SI  NR  for  the  PAMF  is, 


SINR 


4,\*hcpa*  I2 


PAMF  R  C  e 


The  adaptive  pattern  of  the  PA  M  F  i  s, 


PpAMF  I  W  PAMF  e(/ Is’  ftd  ) 


(48) 


(49) 


(50) 


From  the  above  derivation,  we  may  observe  some  features  of  the  PA  M  F  approach: 

•  The  esti  mati  on  of  the  PA  M  F  parameters  (coeffi  ci  ents)  i  s  achi  eved  by  averagi  ng 
measures  not  only  from  a  series  of  range  cells  collected  by  the  same  pulse 
(Equation  (36),  fast-time)  but  also  from  a  series  of  pulses  (Equation  (37),  slow¬ 
time).  FI  encethe  estimation  procedure  is  enhanced  through  the  use  of  both  fast¬ 
time  and  slow-time  data  averaging  process. 

•  It  can  be  seen  from  Figure  1  that  the  PAMF  does  not  use  the  full  covariance 
matrix.  In  another  words,  some  correlations  between  data  collected  by  different 
pulses  are  treated  as  zero.  This  approximation  may  degrade  the  performance  if 
the  actual  correlations  are  high. 

•  The  original  data  snapshot  consists  of  a  time  sequence  of  M  points.  Both  the 
PMF  and  PAMF  processes  are  dimensionality-reduced  transformations  and 
contain  only  a  time  sequence  of  M-p  points.  As  we  know  that  the  Doppler 
resolution  as  well  as  coherent  processing  gain  is  dependent  on  the  length  of  the 
time  sequence,  the  reduction  in  time  sequence  decreases  the  Doppler  resolution 
as  well  as  the  processing  gain.  This  may  not  be  a  concern  if  M  »  p .  FI  owever 
in  the  case  where  the  ratio  ( M-p)/M  is  low,  the  loss  in  the  coherent 
processing  gain  and  the  loss  in  the  Doppler  resolution  may  become  significant. 
Section  7  proposes  a  modified  version  of  PAMF,  which  does  not  induce  any 
dimensionality  loss. 
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5.  Performance  Results 

5.1  Generic  Radar  Model  Data 

A  simple  radar  model  (Ward,  1994)  was  first  used  to  assess  the  performance  of  PMF. 
Table  1  lists  the  parameters  of  the  model.  The  clutter  model,  which  is  often  referred  to 
as  the  first-order  clutter  model,  does  not  take  account  of  temporal  and  spatial 
correlation  effects  (Klemm,  2002)  and  has  the  foil  owing  features, 

•  Constant  cl  utter  coeffi  ci  ent; 

•  U  ncorrel  ated  cl  utter  echoes  from  patch  to  patch; 

•  No  range fol d-over  effect; 

•  No  clutter  intrinsic  motion. 


Table  1:  Parameters  used  in  the  generic  radar  model. 


Parameter 

Specification 

Radar 

Phased  array 

Linear  18-by-4  elements,  element  pattern 
cos(^),  half-wavelength  element  spacing, 
uniform  tapering 

PRF 

300Hz 

N  umber  of  pulses  per  CPI 

18 

Looking  direction 

Broadside,  horizontal 

Platform 

Height 

9km 

Speed 

2va  Kfprfd)  =  1  (DPCA  condition),  where  d 
i  s  the  el  ement  spaci  ng 

Undesired  signals 

Thermal  noise 

Gaussian 

Clutter 

First-order  clutter  model,  CNR  =  47dB 

J  ammi  ng 

Two  equal-intensity  broadband  jamming 
signals  at  azimuth  angles  of  65°  and  130°, 
respectively,  JNR  =38dB 

Under  these  assumptions,  the  covariance  matrix  R;(  can  be  formed  (Ward,  1994,  Dong, 
2005).  Since  R;(  is  strictly  a  Toeplitz-bock-Toeplitz  matrix,  T  =  T(«), 
n  =  0,l,---,M - p-\.  Therefore,  there  is  no  discrimination  between  PMF  and  PAMF, 
and  the  results  of  PMF  and  PAMF  are  the  same.  Figure  2  (a)  shows  the  SI  NR 
comparison  between  theMF  and  PMF  {p  =  3  )  processes.  The  SI  NR  loss  of  PMF  dueto 


13 


DSTO-RR-0313 


dimensionality  loss  compared  to  the  MF  SIN  R  is  shown  in  Figure  2  (b).  The  SIN  R  loss 
duetothedimensionality  lossintheory  is  101og10(15/18)  =  -0.79 dB. 


Figure  2  (a):  5INR  comparison  between  M  F  and  PM  F  (p  =  3).  There  is  a  little  PM  F  51  NR  loss 
compared  to  the  M  F  SIN  R. 


Doppler  frequency  (Hz) 


Figure2  (b):TheSIN  R  lossofPM  F  duetothedimensionality  loss. 

Assume  a  target  signal  \t  which  has  a  constant  amplitude  (SN  R  =  OdB)  and  constant 
Doppler  (lOOFIz)  is  embedded  in  undesired  signals  including  clutter  xc  ~CiV(0,Rc), 
jamming  xy  ~  CA^(0,Ry),  and  thermal  noise  xn  ~  CAT(0,Rn) ,  that  is, 
x  =  xt+xc  +  X  j  +xn  and  R1(  =  R(  +  R  7  +R„ .  Figure  3  shows  the  snapshots  of  clutter, 
jamming,  thermal  noise  and  target  signals.  The  resultant  APMF  is  compared  with 
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AMF  in  Figure  4.  It  can  be  seen  from  the  PAM  result  that  there  is  a  small  loss  in  the 
peak,  and  a  slight  broadening  of  the  mainlobe  due  to  the  dimensionality  reduction 
problem  as  discussed  earlier. 


Figure3:  Snapshots  of  clutter,  jamming,  noise  and  target  signals. 


Doppler  frequency  (Hz) 


Figure 4:  Performancecomparison  between  APMF  and  AMF. 

Finally  the  MF  and  PMF  adaptive  patterns  are  shown  in  Figure  5  assuming  a  target 
signal  to  be  at  zero  azimuth  with  100Hz  Doppler.  The  two  patterns  are  very  similar. 

It  can  be  concluded  from  the  above  analysis,  for  this  simulation,  that  MF  and  PMF 
perform  approxi  mately  the  same  for  thi  s  si  mpl  e  radar  model . 
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Count  of  aamuBi  jngkt 


(b)PMF,  p  =  3 


Figure  5:  M  F  and  PM  F  adaptive  patterns. 


5.2  RLSTAP  Data 

The  second  example  used  for  PAMF  evaluation  was  the  data  generated  by  the  high 
fidelity  airborne  radar  simulation  software,  RLSTAP.  Table  2  and  Table  3  list  radar, 
platform,  envi ronment  and  target  parameters  used  in  the  RLSTAP  simulation. 


Table 2:  parameters  used  in  RLSTAP  simulation. 


Parameter 

Specification 

Radar 

Phased  array 

Linear  20-by-4  elements,  element  pattern 
cos06(^),  azimuth  spacing  0.12m,  elevation 
spacing  0.15m  uniform  tapering  for  transmit 

Carrier  frequency 

1.2GHz 

Polarisation 

VV 

LFM  bandwidth 

2  M  Hz 

PRF 

2  kHz 

N  umber  of  pulses  per  CPI 

32 

Peak  power 

30  kW 

Duty 

10% 

Sample  rate 

0.2  |j,s 

Looking  direction 

30°  off  the  broadside,  horizontal 

Platform 

Height 

7  km 

Speed 

175  m/s 
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Undesired  signals 

Thermal  noise 

Gaussian 

Clutter 

Washington  D.C.  area 

J  ammi  ng 

None 

Table 3:  Target  parameters  used  in  theRLSTAP  simulation. 


Parameters 

Target  1 

Target  2 

T  arget  3 

Height  (km) 

3 

3 

3 

Position  off  main  lobe  direction 

1  km  north 

0.5  km  south 

0 

Radial  velocity  (m/  s) 

-100 

-20 

150 

Doppler  frequency  (Hz) 

-500 

860 

-500 

Range  (km) 

50 

60 

70 

RCS  (sqm) 

1 

1 

1 

RLSTAP  calculates  clutter  returns  based  on  the  United  States  Geological  Survey  Land 
Use  and  Land  Cover  (USGS  LULC)  data  and  the  Digital  Terrain  Elevation  (DTE)  data. 
Details  of  the  clutter  calculation  model  used  in  RLSTAP  are  unknown.  The  clutter 
environment  (Washington  D.C.  area)  seen  by  the  radar  is  shown  in  Figure6. 


Figure  6:  The  LULC  data  superposed  on  the  DTE  data  of  the  1/1/  ashington  D.C.  area. 
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The  unambiguous  range  is  75km  for  the  parameters  given  in  Table  2.  With  the  classical 
4/  3  Earth  radius  model,  the  range  to  the  horizon  is  345km  resulting  in  up  to  four  range 
foldovers.  However  due  to  the  availability  of  the  LCLU  data,  the  simulation  onlytook 
the  first  range  foldover  into  the  account.  That  is,  clutter  returns  from  beyond  150km 
were  ignored.  Because  the  cl  utter  return  decays  as  a  function  of  the  range  to  the  thi  rd 
power,  the  effect  of  higher  range  fold  overs  is  believed  to  be  insignificant.  Accordingly, 
data  collected  by  the  first  pulse  are  discarded  in  the  process,  as  they  do  not  contai  n  the 
clutter  returns  from  ambiguous  ranges.  Therefore  the  actual  number  of  pulses  used  in 
the  process  was  M  =  31. 

To  obtain  the  covariance  matrix  using  the  DL-SM  method,  one  might  think  that  the 
more  range  samples  the  better.  In  fact  this  is  not  necessarily  true  because  samples  are 
generally  not  collected  from  a  statistically  homogenous  clutter  environment.  We  used 
four  different  sampling  methods  and  visually  compared  the  performance  of  MF  for 
target  detection.  The  ranges  bins  used  in  these  four  cases  are, 

•  Case  1:  range  bins  800-2650  with  the  exclusion  of  target  range  bins  themselves 
(the  target  range  bins  for  targets  1,  2  and  3  are  1934,  2266  and  2598, 
respectively)  as  well  as  the  10  nearest  range  bins  on  both  sides  of  each  target 
range  bin; 

•  Case  2:  range  bins  1634-2234  with  the  exclusion  of  1924-1944  (the  target  range 
bin  is  1934) 

•  Case  3:  range  bins  801- 1600  (no  target  in  this  region); 

•  Case 4:  range  bins  1001- 1600  (no  target  in  this  region). 

It  has  been  found  that  case  1  performed  the  worst,  although  it  used  the  most  range 
bins.  Cases  3  and  4  performed  approximately  the  same  and  better  than  the  first  two 
cases.  Case  4  was  used  in  the  following  performance  comparison  between  MF  and 
PAMF. 

The  SI  NR  comparison  between  MF  and  PAMF  {p  =  3 )  is  shown  in  Figure  7  (a).  It  is 
seen  that  the  PAM  F  SI  N  R  shows  some  noticeable  but  not  significant  loss  in  the  clutter 
notch  region  in  comparison  with  the  MF  SI  NR  apart  from  the  loss  caused  by  the 
dimensionality  loss.  The  SI  NR  loss  due  to  the  dimensionality  loss  in  theory  is 

1 0 logI0 (28 / 3 1)  =  -0.44  dB. 
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Figure  7  (a):  SIN  R  comparison  between  M  F  and  PAM  F.  The  PAM  F  SIN  R,  compared  to  the 
M  F  SINR,  shows  some  not  significant  but  noticeable  loss  in  thevicinities  of  the  cl  utter  notch. 


Figure  7  (b):  The  SIN  R  loss  of  PAM  F  due  to  the  dimensionality  loss. 

Results  of  MF  DTSand  PAMF  DTS  are  shown  in  Figure  8.  The  performances  of  the  two 
processors  seem  approximately  the  same  at  first  glance.  Careful  observation,  however, 
indicates  that  the  signal  level  of  target  2  (the  target  in  the  middle  range  bin)  in  the 
PAM  F  process  is  a  few  dB  lower  than  the  level  in  the  M  F  process,  which  is  due  to  the 
PAMF  SINR  I  oss  for  that  frequency  as  shown  in  Figure  7.  Figure9and  Figure  10  replot 
Figure  8  into  2D  plots,  i.e.,  signal  level  versus  range  and  Doppler  frequency, 
respectively.  There  are  only  two  targets  shown  in  Figure  10,  because  target  1  (near 
range)  and  target  3  (far  range)  have  the  same  Doppler  (refer  to  Table  3)  and  are 
superposed.  Figure  9  (a)  shows  the  presence  of  spikes  which  could  increase  the  false 
alarm  rate  or  reduce  the  detection  probability.  We  will  see  in  due  course  that  these 
spikes  in  fact  exist  in  the  results  of  both  MF  and  PAMF  DTS  processes  (though  the 
spikes  are  much  reduced  for  the  PAMF  process  for  this  example  as  shown  in  Figure  9 
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(b))  and  could  become  much  worse  for  MCARM  data.  We  discuss  what  causes  these 
spikes  and  how  to  mitigate  them  in  the  next  subsection. 


‘•CO 


-500 


Ring*  bn 


.1000 


(a)  MF 


Oo«>toftT»q(Ha 


1000 


(b)  PAMF,  p  =  3 


Figure  8:  T  arget  detection  performance  comparison  between  M  F  DT5  and  PAM  F  DTS. 


(a)  MF  (b)PAMF,  p  =  3 

Figure  9:  T  argdc  detection  performance  comparison  between  M  F  DTS  and  PAM  F  DTS,  plotted 
as  signal  level  verses  range  bin. 
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(b)  PAM  F,  p  =  3 


Figure  10:  T  arget  detection  performance  comparison  between  M  F  DTS  and  PAM  F  DTS, 
plotted  as  signal  level  verses  D  oppler  frequency. 


If  we  increase  the  order  of  the  parameter  filter,  p  ,  say,  from  3  to  5,  the  SI  NR  loss  in  the 
vicinity  of  the  clutter  notch  for  the  PAMF  is  reduced  but  at  the  cost  of  Doppler 
resolution  reduction1.  Figure  11  shows  the  detection  result  of  PAMF  with  p  =  5  where 
wecan  see  the  signal  level  of  target  2  having  been  improved  compared  to  the  result  of 
p  =  3  shown  in  Figure 9(b). 


0 


Range  tun 


(a)  Signal  level  versus  range 


Dew**  frequency  (Hz) 


(b)  Signal  level  versus  Doppler 


Figure  11:  U  sing  a  higher  order  parameter  filter,  p  =  5,  improves  the  level  of  target  2  (the 
middle  range  target),  but  at  thecost  of  Doppler  resolution  reduction  in  thePAM  F  processing. 


1  It  means  poorer  Doppler  resolution. 
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5.2.1  Discussions  on  Detection  Test  Statistics 

Detection  test  statistics,  AMF ,  APMF  and  A PAMF  for  MF,  PMF,  and  PAM  F,  respectively, 
have  been  given  in  (9),  (24),  and  (47),  previously.  However  we  have  found  these  test 
statistics  are  sensitive  to  Doppler  frequencies  of  the  mainlobe  clutter,  which  may 
impose  a  difficulty  with  target  detection  and  increase  the  false  alarm  rate.  To  explain, 
let  us  examine  AMF  from  another  perspective.  When  E,t  =1 ,  (7)  can  be  written  as, 


SINRMF^t=\)  = 


I  w  MFe  | 


W  MF 


•  =  e  R“  e 


(51) 


Using  the  result  of  (51)  we  express  AMF  onthedecibel  scaleas, 

101og10  A MF  =  101og10  I  w"Fx  I2  -101og10  SINRMF(%t  =  1)  (52) 


Equation  (52)  indicates  that  on  the  decibel  scale  AMF  is  the  output  of  the  signal  minus 
the  SI  NR  of  the  processor.  In  general,  the  SI  NR  of  the  processor  is  a  constant 
everywhere  except  in  the  frequency  vicinity  of  the  mainlobe  clutter,  where  a  deep 
notch  is  placed  (for  instance,  see  Figure  7).  Therefore  AMF  on  the  decibel  scale  is  the 
output  of  the  processor  plus  a  frequency-independent  constant  (small  in  value)  outside 
of  the  notch  region,  or  a  frequency-dependent  variable  (large  in  value)  in  the  notch 
region.  The  slope  of  the  notch  is  steep  and  varies  rapidly  with  frequency.  On  the  other 
hand,  individual  data  snapshot  x  is  a  random  process  whose  frequency  components  in 
the  notch  region  (i.e.  clutter)  may  vary  significantly.  Therefore,  we  can  expect  that 
fluctuations  of  AMF  in  the  notch  region  could  occur  from  range  bin  to  range  bin,  and 
impose  a  difficulty  on  the  use  of  AMF  for  target  detection.  In  fact  our  interests  in  target 
detection  using  the  MF  process  are  generally  outside  of  the  clutter  notch  region. 
Therefore,  simply  using  the  output  of  the  processor,  101og10 1  w^Fx  |2 ,  as  a  detection 

test  statistic  may  provide  a  better  result.  In  order  to  discriminate,  we  refer  to  this  as  the 
modified  detection  test  statistic  (M  DTS).  Similar  discussion  can  be  made  for  the  other 
two  detection  test  statistics,  APMF  and  Apamf. 

The  detection  results  shown  in  Figure  12  and  Figure  13  are  obtained  with  the  use  of 
MDTS.  It  can  be  seen  that  the  spikes  disappear  while  a  notch  at  the  clutter  Doppler 
frequency  appears. 
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Range  cm 

(a)MFM  DTS,  signal  level  versus  range 
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(b)  M  F  M  DTS,  signal  level  versus  Doppler 


Figure  12:  M  DTS  removes  spikes  in  the  vicinity  of  main  lobe  clutter  D  oppler  frequency,  but 
leaves  a  notch  at  that  frequency. 
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(a)  PAM  F  M  DTS,  p  =  5  ,  signal  level  versus 
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(b)  PAM F  M DTS,  p  =  5  ,  signal  level  versus 
Doppler 


Figure  13:  M  DTS  removes  spikes  in  the  vicinity  of  main  lobe  clutter  Doppler  frequency,  but 
leaves  a  notch  at  that  frequency. 


5.2.2  Minimum  Number  of  Range  Samples 

The  unknown  covariance  matrix  of  undesired  signals  is  usually  estimated  using  range 
samples  as  shown  in  (33).  The  conventional  question  asked  is  how  many  range 
samples  are  required  in  order  to  obtain  a  reliable  estimate  of  the  matrix?  Detailed 
discussions  are  available  in  the  literature  and  are  out  of  the  scope  of  this  report.  Here 
weonly  intend  to  have  some  discussions  from  the  viewpoint  of  practical  application. 
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It  is  well  known  that  given  a  3d  B  tolerance  of  the  optimum  performance  the  number  of 
identically  independently  distributed  (iid)  samples  needed  for  estimating  a  covariance 
matrix  is  twice  the  dimension  of  the  matrix  (Reed  et  al,  1974).  However,  for  the 
equidistant  antenna  array  and  constant  PRF,  the  structure  of  the  covariance  matrix  has 
a  Toeplitz-block-Toeplitz  structure,  and  the  required  minimum  number  of  range 
samples  can  be  greatly  reduced  (Steiner  and  Gerlach,  1998,  Steiner  and  Gerlach,  2000, 
Gerlach  and  Picciolo,  2003,  Bresler,  1998). 

In  general  if  we  know  the  eigenvalue  distribution  of  the  covariance  matrix,  we  can 
estimate  the  minimum  number  of  iid  samples  required  for  the  covariance  matrix.  The 
rank  of  the  covariance  matrix  of  clutter  is  given  by  Brennan's  rule  as  (Brennan  and 
Staudaher,  1992,  Ward,  1994), 

Nr=mt{N  +  j3(M-l)}  (53) 

where  int{-}  denotes  nearest  integer;  j3  =  2va /(fprfd)  is  the  ratio  of  the  normalised 

Doppler  frequency  to  the  normalised  spatial  frequency,  which  is  also  referred  to  as  the 
number  of  the  clutter  fold-overs  (Ward,  1994);  and  d  is  the  azimuth  interval  of  the 
antenna  elements.  In  fact  (53)  does  not  include  the  effects  of  the  range  foldover  and 
clutter  intrinsic  motion.  Therefore,  the  actual  rank  of  the  clutter  covariance  matrix  is 
normally  higher  than  the  theoretical  value  given  in  (53),  and  can  be  determined  by 
examining  the  eigenvalues  of  the  sampled  covariance  matrix.  Figure  14  shows  the 
eigenvalues  of  the  covariance  matrix  Ru  estimated  by  using  range  samples  1001-1600. 
It  can  be  seen  that  the  number  of  significant  eigenvalues  in  this  example  is  about  170. 
Analogous  to  the  general  rule  of  two-fold  iid  samples  of  the  covariance  matrix,  the 
number  of  iid  samples  required  for  the  covariance  matrix  in  this  case  is  then  about 
2x170  =  340. 


Figure  14:  Eigenvalues  of  R„  estimated  using  range  samples  of  1001-1600. 
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One  significant  feature  of  the  PAMF,  according  to  Roman  et  al  (2000)  is  that  the 
method  still  performs  even  using  a  relatively  small  number  of  range  samples  to 

adaptively  estimate  parameters.  This  is  because  the  parameters  A Hf(k), 
k  =  Q,\,- •  •  ,M  —  p ,  are  estimated  using  R;  rather  than  R„ .  The  dimension  of  R;  is 
usually  much  smaller  than  that  of  R„ ,  and  so  is  the  number  of  the  dominant 
eigenvalues.  Therefore  the  number  of  range  samples  required  for  estimating 
parameters  of  PAMF  can  be  significantly  fewer  than  the  number  of  samples  for 
estimating  the  covariance  matrix  in  the  conventional  M  F  method.  Figure  15  shows  the 
eigenvalue  distributions  of  Ry  with  p  =  3  and  p  =  5  estimated  using  range  samples 

1001-1600.  We  make  two  observations  from  the  figure.  First,  the  number  of  dominant 
eigenvalues  are  about  the  same  (26  for  p  =  3  and  28for  p  =  5  )  although  the  dimension 

of  R;  e  €(/,+1),Vx(/,+1)'v  varies  with  p  ,  thus  suggesting  that  using  a  higher  order  of  p 
in  the  PAMF  process  does  not  require  more  iid  samples.  Second,  the  number  of 
dominant  eigenvalues  in  R;  is  significantly  smaller  than  that  of  Ri( ,  so  the  PAMF 
method  requires  far  fewer  samples  in  general. 


Figurel5:  Eigenvalues  of  R^,  p  =  3  and  p  =  5,  estimated  using  range  samples  of  1001- 
1600. 


We  notice  that  the  range  resolution  of  the  system  shown  in  Table  2  is  Ar  =  c/(2B) 
where  c  is  speed  of  light  and  B  the  bandwidth  of  the  LMF.  The  corresponding  range 
sampling  rate  should  be  At  =  1/ B  =  0.5ps  .  The  actual  range  sampling  rate  used  in 
RLSTAP  was  an  over-sampled  rate  of  0.2 ps  .  Therefore,  the  consecutive  range  samples 
may  not  be  iid  samples.  Shown  in  Figure  16  are  the  MF  MDTS  results  with  the 
covariance  matrix  obtained  by  averaging  (a)  300  consecutive  range  samples 
(1301:1600),  and  (b)  200  non-consecutive  range  samples  with  a  interval  of  3 
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(1001:3:1600).  It  can  be  seen  that  the  detection  performance  of  case  (a)  is  starting  to 
deteriorate,  possibly  due  to  only  a  small  portion  of  the  300  samples  being  actually  iid 
sampl  es,  whereas  case  (b)  performs  approxi  mately  the  same  as  before,  even  though  the 
number  of  samples  was  reduced  to  200.  The  sampling  rate  of  case  (b)  is  At  =  3x0. 2 jus  , 
and  the  samples  might  be  considered  as  iid  samples. 


Range  tun 


(a)MFM  DTS  with  300  consecutive  range 
samples 


o 


(00  200  300  400  50)  600  700  800  900 
Range  Bin 


(b)  M  F  M  DTS  with  200  non-consecutive  range 
samples 


Figure  16:  T  arget  detection  results  ofMFMDTS  with  (a)  300  consecutive  range  samples, 
1301:1600  and  (b)  200  non-consecutive  range  samples,  1001:3:1600. 


The  effect  of  using  fewer  samples  in  the  PAM  F  MDTSisshown  in  Figure  17  and  Figure 
18.  In  particular,  the  results  of  Figure  17  and  Figure  18  used  50  (1451:3:1600)  and  30 
(1511:3:1600)  samples,  respectively.  From  Figure  17  we  see  that  for  this  number  of 
samples,  the  MF  process  is  not  functioning,  while  the  PAM  F  process  with  p  =  4  or  5 
still  maintains  a  satisfactory  level,  although  the  PAMF  process  with  p  =  3  starts 
deteriorating.  The  results  of  Figure  18  shows  that  PAM  F  starts  breaking  the  satisfactory 
level  with  p  =  5  ,  whereas  the  result  of  p  =  3  is  unacceptable  because  of  the  detection 
problem  on  the  target  in  the  mid  range. 
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Range  tun 


i 

I 

& 


Range  tun 


(c)  PAMFMDTS,  p  =  4  (d)  PAM  F  M  DTS,  p  =  5 

Figure  17:  Target  detection  results  using  50  (1451:3:1600)  samples:  (a)  M  F  STM  process;  (b)  to 
(d)  PAM  F  M  ST  process  with  p  equal  to 3,  4  and 5,  respectively. 
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Figure  18:  T  arget  detection  using  the  PAM  F  M  DT5  process  with  30  non -consecutive  samples, 
1511:3:160):  (a)  p  =  3  and  (b)  p  =  5. 


5.2.3  Selection  of  Secondary  Data 

Data  used  for  sampling  Ri(  or  R;  are  sometimes  referred  to  as  secondary  data  in 

contrast  to  primary  data  which  areundertest  using  the  statistics  of  the  secondary  data. 
The  primary  data  and  secondary  data  are  mutually  exclusive.  Generally  one  may  use  a 
sliding  window  or  a  fixed  window  to  select  secondary  data  (Ward  and  Kogon,  2004). 
With  the  sliding  window  architecture,  the  CUT  as  well  as  a  few  guard  cells,  which 
normally  sit  in  the  centre  of  the  window,  are  excluded  from  the  sliding  window.  Since 
each  CUT  corresponds  to  a  unique  sliding  window,  it  is  computationally  very 

expensive.  With  the  fixed  window  architecture,  a  single  Ri(  or  R^  estimated  from  a 

fixed  window  is  used  and  applied  to  a  block  of  cells  rather  than  a  single  cell.  All  the 
previous  results  of  RLSTAP  data  were  generated  using  a  fixed  window  as  specified. 
Theoretically,  the  sliding  window  seems  to  be  able  to  provide  better  results  but  has  a 
huge  computational  cost.  However  this  is  not  observed  in  simulation.  Figure  19  and 
Figure  20  compare  the  PAMF  MDST  results  using  a  fixed  window  against  a  sliding 
window.  The  fixed  window  used  was  30  or  50  non-consecutive  samples  with  an 
interval  of  3  (1511:3:1600  or  1451:3:1600).  The  sliding  window  also  used  30  or  50  non- 
consecutive  samples  at  an  interval  of  3,  symmetrically  on  the  sides  of  the  CUT  plus  5 
guard  cells  at  each  side.  That  is,  for  a  CUT  k,  the  sliding  window  consists  of 
(£-6-15x3 :3 : k-6)  and  ( £  +  6:3:£  +  6  +  15x3 )  or  ( £-6-25x3:3:£-6 )  and 
(£  +  6:3:£  +  6  +  25x3).  Weseethat  little  is  gained  using  the  sliding  window  against 
the  fixed  window.  This  is  in  fact  consistent  with  the  findings  reported  previously 
(Dong,  2005).  The  report  (Dong,  2005)  shows  that  first  the  inverse  of  the  covariance 
matrix  is  approximately  invariant  to  changes  in  clutter  range  cells,  and  second,  the 
weighting  vector  is  not  sensitive  to  the  changes  of  grazing  angle  in  a  certain  region. 
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These  two  findings  explain  that  a  single  weighting  vector  can  be  applied  to  a  segment 
of  range  cel  Is. 


PAMF,  p  =  5  ,  with  a  fixed  window  of  30non- 
consecutive  samples 
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consecutive  samples  non-consecuti  ve  samples 


Figure  19:  T arget  d&ection  results  plotted  in  signal-level  versus  range  and  D  oppler  using  the 
PAM  F  ( p  =  5)  M  DTS  process  with  30  and  50  non -consecutive samples. 
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PAM  F  p  =  5  ,  with  a  sliding  window  of  50non- 
conseeuti  ve  samples 


Figure  20:  T  arget  detection  results  plotted  in  signal-level  versus  range  using  the  PAM  F  M  DTS 
( p  =  5)  process  with  30  and  50  non -consecutive  samples. 


5.3  MCARM  Data 

Detailed  descriptions  of  the  multi-channel  airborne  radar  measurement  (MCARM) 
system  can  be  found  elsewhere  (Sloper  etc,  1996,  Fenner  and  Hoover,  1996).  Some  of 
theMCARM  data  analyses  are  also  avail  able  (MITRE,  1999,  RAFDCI,  1999,  Sarkeretal, 
2001).  The  dataset  #5-575  collected  by  the  MCARM  system  was  used  in  this  report.  The 
radar  and  platform  parameters  of  #5-575  are  given  inTable4andTable5,  respectively. 
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T  able  4:  M  CAR  M  radar  parameters 


Frequency/ 

polarisation 

PRF 

CPI 

Pulse 

width 

Duty 

Range 

resolution 

PRI 

(rs/  gates) 

1240  MHz  /VV 

1984 

128 

50.4  ns 

10% 

0.8  |aS 

504/  630 

T able 5:  M  CAR M  platform  parameters 


Height 

Velocity 

Illumination 

Crab  angle 

Antenna  tilt  angle 
from  horizontal  2 

3488  m 

100.1m/  s 

Side-looking 

7.28° 

5.11° 

The  receiver  of  the  MCA  RM  system  consists  of  22  receiving  channels  (modules)  in  two 
azimuthal  rows.  Since  these  two  rows  are  identical,  we  treated  data  received  by  one 
row  as  secondary  data  and  the  data  received  by  the  other  row  as  primary  data.  The 
physical  position  difference  between  the  two  rows  will  in  general  induce  a  phase 
difference  between  signal  arrivals  at  the  two  rows.  This  however  should  impose  no 
effect  as  far  as  generating  the  covariance  matrix  is  concerned.  We  also  only  used  the 
data  collected  by  the  first  50  pulses  in  the  process,  hence,  N  =  1 1  and  M  =  50  in  this 
MCARM  data  study. 

Figure  21  shows  the  clutter  profile  of  the  range  bins  after  compensation  for  the  1  / R3 
range  effect.  The  area  illuminated  by  the  radar  mainlobe  is  mainly  farmland  with 
scattered  houses  (range  bins  200-400  and  500-630,  for  example)  and  bay  water  (range 
bins  400-500,  for  example).  There  is  a  transmitter  signal  leakage  in  range  bin  68  and  the 
useful  clutter  data  are  from  around  range  bin  200  to  the  end. 


2  This  angle  is  the  sum  of  the  antenna  tilt  angle  relative  to  the  platform  plus  the  recorded 
platform  roll  angle. 
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Figure  21:  Clutter  profile  against  range  after  the  1  /  R3  range  effect  is  compensated. 

A  moving  target  at  range  bin  299  has  been  detected  and  reported  previously  (Dong, 
2005).  Range  cells  289  to  309  were,  therefore,  excluded  in  the  secondary  data  for  either 
thefixed  window  or  sliding  window  processing. 

A  fixed  window  consisting  of  range  cells  200  to  630  but  excluding  289  to  309  was  first 
used  to  sample  R„  for  MF  processing  and  for  PAM  F  processing.  The  results  are 

shown  in  Figure  22  to  Figure  24.  Both  processes  are  able  to  detect  the  target  at  range 
bin  299.  Overall,  on  average,  the  noise  floor  of  the  MF  MDTS  result  is  lower.  It  is  also 
worth  noting  that  the  MF  detection  profile  of  signal  versus  range  (Figure  23  (a))  in 
general  follows  the  clutter  profile  (Figure  21)  whereas  the  pattern  of  the  PAMF  MDTS 
(Figure  23  (b))  does  not  show  such  a  feature. 
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Figure  22:  Signal  level  versus  range  and  D  oppler  processed  by  M  F  M  DTS  and  PAM  F  M  DTS 
using  a  fixed  window  of  200:630  with  an  exclusion  of  289:309. 


M  F  M  DTS  PAM  F  MDTS,  p  =  5 


Figure  23:  Signal  level  versus  range  processed  by  M  F  M  DTS  and  PAM  F  M  DTS  using  a  fixed 
window  of  200:630  with  an  exclusion  of  289:309. 
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DocpW  fr^qunrcy  (K?)  Dopplw  IrtOMncy  (H2) 

MFMDTS  PAM F  MDTS,  p  =  5 

Figure  24:  Signal  level  versus  range  and  D  oppler  processed  by  M  F  M  DTS  and  PAM  F  M  DTS 
using  a  fixed  window  of  200:630  with  an  exclusion  of  289:309. 

The  DTS  results  are  shown  in  Figure  25.  We  can  see  from  the  figure  that  both  M  F  and 
PAMF  do  not  sufficiently  suppress  the  clutter,  which  would  significantly  increase  the 
false  alarm  rate  should  the  conventi onal  DTS  be  used. 


Dccpief  (r«qu»ocY  (Hr) 


M  F  M  DTS 


fr*qu*ncy  (Ml) 

PAMF  MDTS,  p  =  5 


Figure  25:  M  F  DTS  and  PAM  F  DTS  processes  produce  many  spikes  at  the  main  lobe  clutter 
Doppler  frequency  deteriorating  the  detection  performance. 


Our  next  intention  is  to  examine  PAMF  with  reduced  secondary  data.  It  is  better  to  use 
iid  samples  for  such  an  intend  on.  The  consecutive  range  cel  Is  of  MCARM  data  may  not 
beiid  samples  as  the  data  were  also  oversampled.  TheLFM  bandwidth  oftheMCARM 
system  was  1MHz,  resulting  in  that  the  range  sample  rate  corresponding  the  range 
resolution  should  be  1.0 jus  whereas  the  actual  range  sample  rate  was  0.8/c  (seeTable 
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4).  Therefore,  non-consecutive  range  samples  with  an  interval  of  2  may  be  considered 
as  iid  samples. 

Apart  from  the  existing  target  in  range  bin  299,  we  also  injected  a  target  signal  into 
range  bin  500.  The  target  has  a  constant  amplitude  of  30dB  below  the  mean  clutter  of 
that  bin  and  a  constant  Doppler  frequency  of  -200Hz.  The  results  of  MF  MDTS  and 
PAMF  MDTS  using  three  fixed  windows,  200:2:600,  200:2:300  and  200:2:260  (excluding 
the  block  of  289:309  if  applicable)  are  shown  in  Figure  26,  Figure  27  and  Figure  28, 
respectively.  It  can  be  seen  that  the  PAMF  process  is  more  robust  in  terms  of  using 
reduced  sample  data.  For  the  M  F  process,  the  injected  target,  which  is  detectable  with 
the  use  of  a  large  number  of  sample  data,  becomes  non-detectable  when  the  number  of 
samples  is  reduced.  For  the  PAMF  process,  on  the  other  hand,  the  injected  target  is 
detectable  with  all  three  fixed  windows  used.  The  original  target  in  range  bin  299 
however  is  always  detectable  even  for  the  MF  process  with  the  window  of  200:2:260. 


Ranga  On  Rang?  on 


MFMDTS  PAMF  MDTS,  p  =  5 

Figure  26:  T  arget  detection  using  MF  MDTS  and  PAMFMDTS  processes  with  a  fixed 
window  of  200:2:600  (exclusion  of  389:309). 
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Figure  27:  T  arget  detection  using  M  F  M  DTS  and  PAM  F  M  DTS  processes  with  a  fixed 
window  of  200:2:300  (exclusion  of  289:300). 


Range  cm  Range  on 


MFMDTS  PAM F  MDTS,  p  =  5 

Figure  28:  T  arget  detection  using  M  F  M  DTS  and  PAM  F  M  DTS  processes  with  a  fixed 
window  of  200:2:260. 

Finally  the  results  of  PAM  F  M  DTS,  with  the  use  of  sliding  windows  of  50  and  30  non- 
consecutive  samples,  are  shown  in  Figure  29.  It  shows  again  that  compared  to  the 
results  of  the  fixed  windows  shown  in  Figure  27  and  Figure  28,  there  is  no 
improvement  in  terms  of  target  detection,  although  a  significant  amount  of  additional 
computation  is  needed  for  implementation  of  the  sliding  window. 
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Rang*  Cun 

(a)  PAMF  MDTS,  p  =  5  ,  with  a  sliding 
window  of  50  non-consecutive samples 


S 

I 


Rang* bm 

(b)  PAMF  MDTS,  p  =  5  ,  with  a  sliding 
window  of  30  non-consecutive  samples 


Figure  29:Target  detection  using  the  PAM  F  M  DTS  process  with  a  sliding  window  of  (a)  50 
non-consecutive  samples  and  (b)  30  non-consecutive  samples. 


5.4  Conclusions 

The  performance  of  the  PAMF  approach  has  been  assessed  in  this  Section  using  three 
types  of  radar  data.  The  results  of  the  conventional  MF  method  have  served  as 
benchmarks  in  the  assessment.  Based  on  the  results,  we  have  observed  the  following 
points: 

•  MF  outperforms  PA  M  F  if  sample  data  are  ample. 

•  PAMF  is  much  more  robust  if  sample  data  are  limited. 

•  The  sliding  window  processing  does  not  seem  to  show  a  better  detection 
performance,  and  the  fixed  window  processing  is  satisfactory  and  very 
computationally  efficient. 

•  The  conventional  DTS  by  the  inclusion  of  numerically  unstable  data  in  the 
clutter  notch  can  produce  false  alarms  and  a  modified  version  (MDTS) 
generally  outperforms  the  conventional  DTS. 

Operational  counts  for  the  PAMF  algorithm  and  its  computational  savings  compared 
to  the  conventional  MF  algorithm  are  discussed  in  the  next  Section. 


6.  Estimation  of  Operational  Counts 


This  section  estimates  operational  counts  (ops)  for  the  MF  and  PAMF  algorithms.  In 
the  estimation,  we  only  count  complex  multiplication/ division,  and  ignore 
addition/ subtraction.  For  simplicity  we  also  assume  that  no  special  digital  signal 
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processing  (DSP)  hardware  is  used,  so  that  all  calculations  are  treated  with  the  same 
weight.  Special  structures  of  matrices  such  as  complex  conjugate,  which  often  lead  to 
special  algorithms  to  save  ops,  are  also  not  considered.  For  instance,  the  number  of  ops 

for  a  matrix  multiplication  A  =  BC  is  Imn  where  A  e  <C'"X" ,  B  e  <C'"x/  and  C  e  <C/x" . 

Ops  required  for  the  MF  and  PAMF  algorithms  are  listed  in  Table  6  and  Table  7, 
respectively. 


T able  6:  0  ps  required  for  the  M  F  process. 


Algorithm 

Task 

Reference 

Eqn 

Ops 

CMF{  1) 

K  H 

R„  = 

k= 1 

(36) 

K(MN)2 

CMf(  2) 

(1) 

c MAO3 

CMf(  3) 

Wmf  =R«'e 

(1) 

(. MN )2 

CMf(  4) 

^MF  0  =  WMFe 

(9) 

MN 

Cmf  (5) 

Iw^  xl2 

A  I  VVMFA  1 

(9) 

MN 

1VMF  _  . 

AMF  0 

N  ote  that  the  actual  ops  for  CMF  (1)  may  be  reduced  by  a  half,  if  we  take  account  of  the 
Hermitian  structure  of  R(( .  Similarly,  the  ops  required  for  CMF(  2)  may  be  reduced  by 

a  half  as  R^1  is  also  Flermitian.  The  ops  required  for  the  inverse  of  a  matrix  arejustified 
by  Isaacson  and  Keller  (1966). 
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T  able  7:  0  ps  required  for  the  PAM  F  process. 


Algorithm 

Task 

Reference 

Eqn 

Ops 

C PAMF  (1) 

kf 

(35)  and 
(36) 

KN2 [M 2  -  (M  -  p)2  +(M  -  p)] 

CpAMF  (2) 

Af 

(32) 

( pN )3  +  p2N 3 

CpAMF  (3) 

T 

(37)  and 
(39) 

[(/>  +  l)2+(/>  +  l)  +  l]/V3 

C PAMF  (4) 

u 

(40) 

pN 2  ( M  -  p) 

C PAMF  (5) 

s 

(40) 

N2(M  -  p) 

C PAMF  (6) 

£ 

(40) 

pN2  (M  -  p) 

C PAMF  ("7) 

V 

(40) 

N2(M  -  p) 

CpAMF  (8) 

A PAMFO  =  S  S 

(40) 

N(M-p ) 

C PAMF  (9) 

.  1  s^v  I2 

^ PAMF  -  . 

A  PAMFO 

(40) 

N(M-p ) 

In  the  MF  process  assuming  the  Doppler  resolution  used  in  Algorithm  CMf( 3)  is 
1/(2 M)  (normally  the  number  of  Doppler  bins  is  equal  to  M ),  we  need  2 M 
repetitions  of  CMF( 3),  CMF( 4)  and  CMF(5)  in  searching  for  possible  target  signals 
with  unknown  Doppler  frequencies  and  amplitudes  in  a  single  range  cell.  If  we  use  a 
fixed  window  for  sampling  RM  the  total  ops  for  processing  a  segment  of  Klg  range 
cel  Is  collected  by  a  CPI  are, 

Cmf  =  Cmf  (\)  +  Cmf(2)  +  2  M\Cmf  (3)  +  CMF  (4)]  +  2  KrgMCMF  (5)  (54) 

Similarly  assuming  that  the  Doppler  resolution  used  in  calculation  for  the  PAMF 
processing  is  1/(2M),  and  a  fixed  window  is  used  for  sampling  R  f -,  then  the  total  ops 
for  processi  ng  a  segment  of  Krg  range  cel  I  s  col  I  ected  by  a  C  PI  are, 


C 


PAMF 


-  C PAMF  (1)  +  CpAMF  (2)  +  CpAMF  (3)  +  2M\Cpamf  (4)  + 

CpAMF  (5)  + 

CpAMF  (8)]  + 

Krg  [CpAMF  (6)  +  CpAMF  (7)]  +  2KrgMC  PAMF  (9) 


(55) 


Shown  in  Figure  30  is  the  ratio  of  the  PAMF  ops  to  the  MF  ops,  which  is  a  function  of 
the  number  of  pulses  M  ,  the  number  of  antenna  channels  N  and  the  number  of  range 
cells  Krg .  The  other  two  parameters  used  in  the  calculation  include  K  =  50  and  p  =  5  . 
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Thefixed  window  processing  isassumed.  Itcan  be  seen  that  the  computational  savings 
of  the  PAMF  process  varies  with  parameters.  For  a  case  where  M  and  N  are  large, 
and  K  small,  the  computational  advantage  of  the  PAMF  process  is  significant.  On 

the  other  hand,  if  M  and  N  are  small,  and  K  large,  there  is  not  much 
computational  advantage i n  the  PAMF  process.  This  is  becausewhen  Krg  is  small,  the 
ops  for  the  inversion  of  Ri(  is  dominant  in  the  MF  processing.  However  with  an 
increase  in  Krg ,  the  portion  for  computing  the  inverse  of  R„  in  the  whole  process 
becomes  I  ess  domi  nant. 
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(a)  Krg  =  100 


(b)  Krg  =  500 


Figure30:  Ratio  of  thePAM  F  ops  to  theM  F  ops,  asa  function  ofthenumber  of  pulses  M ,  the 
number  of  antenna  channels  N  and  thenumber  ofrangecells  K  . 


Ops  required  by  the  MF  and  PAMF  processes  are  shown  in  Figure  31.  It  shows  the 
amount  of  ops  even  for  the  PAM  F  process  is  still  huge,  considering  the  time  frame  for  a 
CPI  to  bein  a  scale  of  milliseconds  for  airborne  radars.  Therefore  how  to  further  reduce 
ops  will  still  bean  active  research  area  in  the  future.  The  figure  also  shows  that  for  the 
same  process,  an  increase  in  the  number  of  range  cells  to  be  processed  does  not  incur  a 
significant  increase  of  ops.  Therefore,  using  a  fixed  window  to  process  a  large  chunk  of 
range  cel  I  s  shoul  d  be  a  preferabl  e  opti  on. 
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- MF,  N  =  10 

MF,  N  =  30 
PAMF,  N  =  10 
PAMF,  N  =  30 


(a)  Krg  =100 


(b)  Krg  =500 

Figure 31:  Ops  required  for  theM  F  and  PAM  F  processes. 


7.  PAM  F  with  Combined  Forward  and  Backward 

Coefficients 

One  of  the  problems  of  the  PAMF  method  is  the  dimensionality  loss  mentioned  in 
Section  3.  Using  the  pi\v  order  AR  parameters  A Hf(k),  k  =  \,---,p,  we  can  only 
estimate  [x(n)\n  =  p,p  +  l,---,M -\}  from  the  original  data  sequence 
{x(/?)| n  =  0, 1, •  •  • , M  —  l} .  It  means  that  the  original  M  -point  time  sequence  now 
becomes  an  (M  -  p)  -point  time  sequence.  As  we  know  that  both  the  Doppler 
resolution  and  processing  gain  are  proportional  to  the  number  of  points  in  a  time 


42 


DSTO-RR-0313 


sequence,  the  dimensionality  reduction  directly  results  in  a  loss  in  the  Doppler 
resolution  and  processing  gain.  The  impact  can  be  small  or  significant  depending  on 
the  values  of  M  and  p.  In  the  case  where  the  ratio  of  (M-p)/M  is  small,  the  impact 
can  be  significant. 

Equation  (14)  is  a  process  of  the  forward  linear  prediction  which  estimates  the 
sequence  {x(n)|n  =  p,p  +  l,---,M -l}.  The  sequence  {x(n)|«  =  0,1, which 

cannot  be  estimated  using  the  forward  linear  prediction,  however,  may  be  estimated 
with  the  use  of  the  backward  linear  prediction  as  (Marple,  Jr,  1987), 

x(n)  =  -ZAf  (k)x(n  +  k)  n  =  0,1, --^p-l  (56) 

k=l 

where  A f(£),  k  =  l,--,p,  are  backward  prediction  parameters.  Together  with  the 

forward  prediction  given  in  (14),  wenow  haveafull  dimension  of  thelinear  prediction 
residual  as, 


s(n)  =  x(n)  -  x(n) 


ZAf  (k)x(n  +  k)  n  =  0, 1,  -  —  1 

k=0 

X  A Hf(k)x(n-k)  n  =  p,p  +  \,---,M -1 
j=o 


(57) 


The  remaining  process  is  the  same  and  there  is  no  need  to  repeat  it  here  except  to 
define  the  estimates  of  the  backward  parameters,  A b(k) ,  k  =  1, •••,/>.  In  order  to  find 
these  parameters,  we  need  to  find  the  mini  mum  squared  residual  given  by  (58), 


H  P~l  H  M~X  H 

min{£  £}  =  minX£  (»)£(«)  + min  X  £  («)e(«) 


n= 0 


n=p 


(58) 


Thefunction  min^E}  in  (58)  isexpressed  as  the  sum  of  two  terms  indicating  that  the 
backward  parameters  and  forward  parameters  are  uncorrelated  and  can  be  separately 
determined.  The  estimate  of  the  forward  parameters  given  in  (32)  therefore  remains 
unchanged.  The  estimate  of  the  backward  parameters  is  given  below  as  (following  a 
similar  mathematical  manipulation), 


A b  ~ 


A,(0) 

A,(l) 

A  h(p) 


LN 


(59) 
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where  Rhh  eCpNxpN  and  Rft0  eCpNxN  are  sub-block  matrices  of  R6  e 
given  as, 


R.= 


o 

o 

Pi 

1 _ 

o 

-C) 

Pi 

_ 1 

I  R  „(k:k  +  p) 


k= 0 


(60) 


It  is  worth  noting  that  the  order  of  Ah  in  (59)  is  from  low  to  high,  contrary  to  the  order 
of  Af  in  (32)  which  is  from  high  to  low.  Also  the  divisions  of  the  sub-block  matrices  in 
(33)  and  (60)  used  for  computing  Af  and  Ah ,  respectively,  aredifferent. 


The  manipulation  matrix  B  given  in  (20)  should  be  modified  accordingly  to  include 
the  backward  prediction  as, 


IjV  A b  ( P ) 

1 N  A b  ( P ) 

IjV  A b  ( P ) 


A-f(p)  •••  Ijv 

A  Hf(j>)  -  lN 

A  Hf(p)  -  IN 


(61) 


Similarly  the  corresponding  mean  temporal  residual  covariance  matrix  is  given  as, 


R 


£ 


)R 


bs 


v  M-p 


i  R 


fe 


where 


(62) 


Rto=lzRt(«)  =  lAfRtAt 

P  n= 0  p 


(63) 


Aimed  at  avoiding  the  dimensionality  loss  problem  we  have  modified  the  PAMF 
algorithm  using  both  the  forward  and  backward  predictions  to  maintain  the 
dimensionality  of  the  time  sequence.  The  above  formulation  uses  p  points  of 
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backward  prediction  and  M-p  points  of  forward  prediction.  Obviously  we  can  form 
other  combinations  if  desired.  We  could  select  a  combination  containing  a  maximum 
number  of  points  including  M-p  points  of  backward  prediction  and  M-p  points  of 
forward  prediction,  for  instance.  However,  according  to  the  linear  algebra,  the 
maximum  number  of  independent  vectors  after  a  linear  transform  is  no  greater  than 
the  original  number  of  the  independent  vectors.  Therefore,  we  do  not  expect  that  there 
woul d  exi  st  a  combi  nati  on  whi ch  contai  ns  more  i  ndependent  vectors  than  the  ori  gi  nal . 

To  view  the  impact  on  the  Doppler  resolution  and  the  processing  gain  due  to  the 
dimensionality  reduction  of  PAMF,  the  example  shown  in  Figure  4  is  re-examined.  To 
clearly  show  the  impact  we  have  chosen  M  =  10  and  p  =  5  .  The  results  are  shown  in 
Figure  32.  The  Doppler  resolution  of  PAMF  with  forward  prediction  for  this  example  is 
only  a  half  of  the  Doppler  resolution  of  the  conventional  STAP,  since  the  PAMF 
process  loses  50%  of  the  time  sequence.  On  the  other  hand  we  seethe  resolution  of  the 
PAM  F  with  both  backward  and  forward  predictions  is  the  same  as  that  of  the  STAP, 
because  of  no  dimensionality  loss  of  the  algorithm.  A  3dB  loss  in  the  signal  peak  level 
for  the  process  of  the  PAM  F  with  forward  prediction  only  is  also  observed. 


Cbppler  frequency  (Hz) 


Figure  32:  Doppler  resolution  comparison  between  thePAM  F  with  forward  prediction  and  the 
PAM  F  with  both  backward  and  forward  predictions.  The  ST A  P  result  is  also  shown  as  a 
benchmark. 

As  the  ratio  of  (M  - p)/M  increases,  the  impact  of  the  dimensionality  loss  fades.  The 
detection  results  for  target  1  in  the  RLSTAP  dataset  is  plotted  in  Figure  33.  In  this 
example  the  values  of  M  =  31  and  p  =  5  result  in  a  16%  decrease  in  the  Doppler 
resolution  if  the  PAMF  with  forward  prediction  is  used. 
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Doppler  frequency  (Hz) 


Figure33:  Theimpact  of  the  dimensionality  loss  on  theDoppler  resolution  and  processing  gain 
fades  as  theratio  of  ( M  -  p)/M  increases.  Shown  is  thetarget  detection  of  target  1  in  the 
RLSTAP  dataset  with  M  =31  and  p  =  5. 


8.  Summary 

M  uch  research  effort,  notably  that  carried  out  and  sponsored  by  AFRL,  has  been  made 
in  developing  and  evaluating  the  PAMF  approach  for  the  last  decade.  This  report  has 
traced,  analysed  and  evaluated  this  technique  in  great  detail.  Accessible  publications  of 
PA  M  F  often  omit  detai  I  s  of  the  mathemati  cs  due  to  space  I  i  mi  tati  on,  maki  ng  the  papers 
difficult  to  follow.  This  report  therefore  has  independently  provided  full  mathematical 
derivation  for  a  full  understanding  of  the  algorithm.  The  performance  of  the  approach 
has  then  been  assessed  using  three  different  types  of  data,  namely  data  generated 
using  a  simple  generic  radar  model,  data  generated  by  high  fidelity  radar  simulation 
software,  RLSTAP,  and  real  radar  data  collected  by  the  MCARM  system.  In  the 
assessment,  the  results  of  the  conventional  space-time  MF  process  have  been  used  as 
benchmarks  for  comparison.  The  ops  for  the  approach  have  also  been  estimated  and 
computational  savings  of  PAMF  in  comparison  to  the  MF  algorithm  have  been  given. 
Finally  a  modification  to  the  PAMF  approach  is  proposed.  The  modified  PAMF  uses 
both  the  forward  and  backward  prediction  parameters  to  maintain  the  dimensionality 
of  the  process,  so  that  the  Doppler  resolution  and  processing  gain  is  maintained. 

Based  on  the  study,  the  PAMF  algorithm  has  two  advantages  over  the  conventional 
MF  algorithm: 

•  The  ops  requirement  is  less.  Computational  savings,  compared  to  the  MF 
algorithm,  vary  from  less  than  20%  for  a  system  with  a  small  number  of  pulses 
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M  and  a  small  number  of  channels  N ,  and  a  large  number  of  range  bins  Krg , 
to  more  than  90%  for  a  system  with  large  M  and  N ,  and  small  Krg . 

•  The  PAMF  algorithm  is  much  more  robust  than  the  MF  algorithm  in  the  case 
w  here  the  set  of  samp  I  e  data  i  s  smal  I . 

Problems  of  the  PAMF  process  include 

•  If  PAM  F  uses  the  forward  prediction  parameters  only,  the  equivalent  temporal 
independent  measurements  are  reduced  from  M  to  M  - p  ( p  is  the  order  of 
the  filter,  and  M  is  the  number  of  pulses  in  a  CPI),  which  in  turn  reduces  the 
Doppler  resolution  and  the  processing  gain.  The  impact  can  be  significant  if  the 
ratio  of  (M-p)/M  is  small.  FI owever,  as  proposed,  if  PAMF  is  formed  using 
both  the  forward  and  backward  prediction,  thereis  not  a  loss  in  dimensionality. 

•  The  output  of  the  PAM  F  M  DTS  processing  does  not  seem  to  follow  the  pattern 
of  the  original  clutter  profile  whereas  the  MF  MDTS  does  (see  Figure  23  and 
Figure  21).  Whether  this  imposes  a  problem  to  target  detection  warrants  a 
further  investigation. 

A 1 1  si  mul  ati  ons  have  suggested  that  usi  ng  the  si  i  di  ng  wi  dow  architecture  to  sampl  e  the 
covariance  matrix  for  each  range  bin  does  not  statistically  improve  results.  Therefore  a 
preferable  option  is  to  use  the  fixed  window  architecture  to  sample  the  covariance 
matrix  and  process  a  block  of  range  bins,  which  significantly  saves  computational  cost. 

Overall  the  amount  of  ops  required  for  the  PAMF  algorithm  is  still  huge  considering 
the  time  frame  available  for  airborne  radars.  Flow  to  further  reduce  ops  demand  in 
processing  phased  array  data  warrants  future  research. 
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